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ABSTRACT 


The present study has been undertaken to determine 
numerically the flow structure and heat transfer characteristics 
in a rectangular channel with built— in vortex generators in the 
form of slender delta wing and winglet pair. Each wing or winglet 
pair induces streamwise longitudinal vortices behind them which 
cause a spiralling flow and entrain fluid from outside into the 
core of the vortices. Finally, these vortices disrupt the growth 
of thermal boundary layer and bring about the enhancement of heat 
transfer between the flowing fluid and the channel walls. The 
prediction of complex three dimensional flow structure and 
augmentation of heat transfer are made by making use of the 
solution of complete Nav i er-Stokes and energy equations. The 
geometrical configuration represents either an element of a 
compact gas-liquid fin-tube crossflow heat exchanger or an element 
of a plate-fin crossflow heat exchanger. ' 

In practice the vortex generators can be mounted on the 
flat surfaces of the above mentioned heat exchangers by punching 
or embossing the flat surfaces. The effect of these punched holes, 
known as stamping, on heat transfer and skin friction 
characteristics has also been accounted. The effect of angle of 
attack of the wing generators, their geometrical shape and 



Reynolds number on heat transfer enhancement and skin friction 

I 

distribution have been obtained as well. It is found that the 
enhancement in heat transfer is more for the case with built-in 
vortex generators and without any stamping than for the case with 
punched holes beneath the vortex generators. Heat transfer 
coefficient is found to increase with increase in angle of attack. 
Finally, the delta wing shows a better performance in heat 
transfer enhancement as compared to winglei pair for the same 
geometrical area of the wing and the winglet pair. 

However, most convective heat transfer processes 
encounter two types of losses, namely, losses due to fluid 
friction and those due to heat transfer across finite temperature 
gradient. These two phenomena are manifestation of thermodynamic 
irreversibility and an evaluation of the process from this 
standpoint is referred to as the second-law analysis. A second^law 
analysis is made for the above mentioned augmentation techniques 
and conclusions are drawn about the influence of type of vortex 
generators < wings / winglets ) on irreversibility. 
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NOMENCLATURE i 

B = channel width 

Br = ratio of wing span to width of the channeHb/B) 

b = wing span 

= skin friction, [z C I7 )J/ ( 

= combined spanwise average friction coefficient 
equation (2?) 

Cp = specific heat of the fluid 

Ec = Eckert number 

i 

H = channel height 

h = heat transfer coefficient, -kf -3— "1 /<T -T. ) 

1. ffy ■'v/ V b 

I = rate of destruction of exergy in the channel 

k = thermal conductivity of the fluid 

M = merit function 

Nu = local Nusselt number based on bulk 

temperature of the fluid 
Ru = spanwise average Nusselt number 

TJu = average Nusselt number for the entire channel 

e 

Ns = total rate of nond imens 1 onal entropy 

generation in the entire channel , equation <1?) 



= rate of nondimens i onal entropy generation per 
unit volume, equation (12) 

= nondimens i onal pressure, p/ pU 

'CIV 

= pressure 

= Prandtl number, pC / k 

P 

= total rate of heat transfer to the fluid 
= total rate of exergy transfer corresponding to 
energy transfer Q 
= wall heat flux 
= Reynolds number, U H / v 

av ' 

= ratio of wall temperature to ambient 
temperature =T /T -T /T 

V ' a w ' CO 

= total rate of entropy generation in the entire 
channe 1 


= rate of entropy generation per unit volume, 
equation <8> 

= temperature / 

= ambient temperature 

= average bulk temperature for the entire channel 
= time 


= u/ U 

' O-V 

= axial component of velocity 


V 
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V 

W 

w 

X 

X 

Y 

y 

z 
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= vertical component of velocity 


= spanwise component of velocity 
= x/ H 

= axial dimension of coordinates 
= y/ H 

= vertical or normal dimension of coordinates 

= z/ H 

= spanwise dimension of coordinates 


reek symbols 

« = aspect ratio of the channel, B / H 

A = divergence of velocity 

6.. = Kronecker delta 

/J = dynamic viscosity of the fluid 

V = kinematic viscosity of the fluid 

O' = stress tensor 

9^ = average nondimens ional bulk temperature for 

the entire channel = <Tk-T )/(T - T> 

T = nondimensional time 

3 ? # 

A = aspect ratio of the wing, b'/ S 



Subscripts 


w = wa 1 1 

b = bulk condition 

av = average 

1 = bottom fin-plate 

2 = top fin-plate 

sa = spanwise combination of top and bottom plate 
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CHAPTER I 


INTRODUCTIOrJ 


Longitudinal vortices embedded in shear flows occur in 
many flow situations and have remarkable influence in design and 
development of heat exchangers. Heat exchangers are used 
extensively in the energy industry, the petrochemical industry, 
the process industry and the food industry. Improvement in heat 
exchanger design would have a direct impact on material saving, 
space saving and efficient use of energy. Figures l<a> and l(b> 
show typical arrangements of two different heat exchanger cores. 
It has been seen that the heat transfer coefficients on the flat 
surfaces of the above mentioned heat exchangers are significantly 
low. Augmentation of heat transfer is of special interest in 
these types of channel flows where the heat transfer between the 
fluid and the channel walls deteriorates as the boundary layer 
grows on the channel walls and the flow tends to become fully 
developed. 

In order to disrupt the growth of thermal boundary layer 

I 

and, thereby, increase heat transfer between the gas and the fin 
in the case of gas-liquid fin-tube crossflow heat exchanger and 
between the flowing fluid and the plates in the case of plate-fin 



Liquid 



(a) 


Fig 1 Typical arrangement of heat exchanger cores 

(a) Gas- liquid fin-tube crossflow 

(b) Plate-fin (single or multi pass') 



Fig. 2 Protrusions in the form of (a) delta wing 
and (b) winglet pair on the flat surface 
to enhance heat transfer. 


heat exchangers Protrusions could bp ^^^un^pd on the flat surfaces. 
The protrusions in the form of slender wing~type vortex generators 
C Fig. 2 <a) and (b) ] can bring about the desired augmentation of 
heat transfer at the expense of relatively less increase in 
pressure drop. These wing generators will take fluid from the 
underside of the wake and swirl it around the upper side 
entraining fluid from the channe 1 -wa 1 1 -regi on to the 
vortex-core-^regi on . This mechanism finally culminates in the 
disruption of the growth of thermal boundary layer and consequent 
increase in heat transfer coefficient. Because of favourable 
pressure gradient in the channel, the longitudinal vortices will 
be stable and their effect will be perceived over an area which is 
many times the size of the wing generators. From the consideration 
of manufacturing, these vortex generators can easily be fabricated 
by cutting delta shape on the plates and then pushing them up. 
Finally, the slender wing generators can also act as spacers for 
the plate fins. Typically, the Reynolds numbers for such 
applications lie below 5000. Experimental investigations by 
Edwards and Alker Cl], Russel, Jones and Lee [2] and Fiebig et al 
C5] and computational analysis by Fiebig al [43 and Biswas et 
al C?3 show that the vortex generators in the form of slender 
delta wings and winglet pairs are indeed very effective in 
augmentation of heat transfer. Velocity field and heat tr4nsfer 
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measurements were reported for slender wings embedded in a 
turbulent boundary layer <see Eibecl' and Ea^on [63; Ligrani et al 
£7J >. Recently, a computational study has been performed by 
Biswas and Chattopadhyay C8] for lammnr flows in a geometrical 
configuration of delta wing placed inside a channel. They have 
discussed about the influence of angle of attack and Reynolds 
number on heat transfer and skin friction coefficients. Fiebfg et 
al C93 have shown that the delta wing^ are more effective than 
the winglet pairs with regard to augmentation of heat transfer 
coefficient on a flat surface. It has also been shown that the 
drag generated in the channel due to winglet pair is considerably 
less than that due to the wing. In this investigation we shall try 
to compare the performance of wings and winglets by taking into 
consideration both the parameters, nannoiy, the augmentation of 
heat transfer and the associated flovj losses. 

The aforesaid comparison may be accomplished with the 

i 

help of second-law analysis. The inherent irreversibility 
associated with any convective heat transfer process results in 
penalty to useful power. Estimation of irreversibilities of a 
process from the entropy generation rate and the minimization of 
the same by a suitable adjustment of the governing input 
parameters is usually referred to as second-law or thermodynamic 
optimization of the process. In such analyses, the 
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irreversibilities due to external interaction of energy and 
internal dissipative effects are simultaneously taken care of. 

The present work provides a qualitative and quantitative 
performance data for delta wing and winglet pair in a channel for 
different operating conditions. Solution of complete Nav i er-Stokes 
and energy equations and a detailed analysis of flow structure 
provide deep insight into the mechanistic modelling of the 
complicated flow configuration. Finally, a numerical evaluation of 
entropy generation in the heat transfer process for the proposed 
model is presented. Based on an evaluation of irreversibility, 
this study draws some conclusions about the influence of type of 
vortex generators < wings / winglets ) on efficient use of energy. 
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t 


t 



CHAPTER II 


LITERATURE SURVE Y 

The objectives of the present study has been briefly 
outlined in the previous chapter. The prf»sont chapter deals with a 
survey of relevant literature in the field o f> augmentat i on of heat 
transfer by surface mounted obstacles. The survey helps in 
suggesting further works that should be carried out to accomplish 
the objectives enumerated earlier. 

In both gas-liquid fin-tube crossflow heat-exchangers 
and plate-fin heat exchangers (single or multi-pass), the heat 
transfer coefficient between the flowing fluid and the flat plates 
are very low. In order to improve the heat transfer performance of 
the heat exchangers, protrusions could be mounted on the flat 
surfaces. Of course, the heat transfer enhancement with the help 
of obstacles is always accompanied by a corresponding rise in 
pressure drop. It has been observed, however, that the protrusions 
in the form of slender delta wing or winglet pair can augment heat 
transfer significantly at the expense of relatively less increase 
in pressure penalty. The enhancement is obtained due to the 
generation of streamwise longitudinal vortices behind the 
protrusions and consequent better mixing of hot and cold fluid in 
the downstream. The slender wings and winglets are streamline 


I 



bodies, hence the associated pressure drop is also less. Van Dyke 
CIOD has illustrated the flow v i sus 1 i Jta 1 1 on results of the 
vortices formed by delta wing and slenHor bodies. One of the flow 
visualization is reproduced in fig.2fc>. 

The first experimental result leqnrdmg augmentation of 
heat transfer by vortex generators has been reported by Edwards 
and Alkar C13. They, however, used bluff bodies (cubes) as vortex 
generators. Russel, Jones and Lee C2] has presented some 
quantitative results for heat transfer and pressure loss due to 
small aspect ratio winglets. They have compared square winglets 

I 

with delta wings and found, in contrary to the results of other 
investigators, a better performance for the winglets. 

Experimental investigations due to Fiebig, Kallweit, and 
Mitra C33 is the first systematic studv to compare the performance 
of different kinds of vortex generators m the reynolds number 
range of 1560 to 2270. Their observations have established the 
fact that the delta wing is the best vnrtev generator from heat 
transfer point of view. Another important feature of their 
observations is that the heat transfer roefficient increases with 
increase in angle of attack. They have also pointed out that the 
ratio of highest heat transfer enhancement to the friction factor 
increment is achieved with delta wing at a small angle of attack, 
about 10*. Their flow visualization by' Laser technique has 
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Fig. ace) Plan view of vortices formed by a slender delta wing. 

couitesy of II. Werlc, ONIiRA, I'lancc 
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revealed that the concentrated vortex rair generated by a small 
aspect ratio delta wing at large angle of attack has elliptic 
shape due to the influence of channel walls. High aspect ratio 
delta wings have shown a turbulent struc^'ure m the vortex core. 

Formation of streamwise longitudinal vortices behind a 
slender aerodynamic object is a research topic of considerable 
interest for several years. Both the theoretical and experimental 
investigations on flow past a delta wing have been conducted and 
reported in literature by a number of researchers. Among them 
Winter Cllil, Fink C123, Mardsen et al [153, Peckham C143. Earnshaw 
et al C15] and Lawford C163 are noteworthy. 

Experimental measurements by Winter was followed by the 
very early nonlinear theory of rectangular wings of Bollay C173. 
More realistic representation of the physics of vortex formation 
has been furnished by Brown et al [183, Mangier et al [193 and 
smith et al [203. Use of thin aerofoil theory provides good 
approximation of the the strength of the leading edge vortices and 
the lift. Of course, in trailing edge part, their results show 
some discrepancy with experiments. Rehbach [213 and Kandil et al 
[223 have obtained better results with non-conical theory. 
Computational solution of non-conical flow field and pressure 
distribution on wing surface has also been attempted by Rubbert 
and co-workers [25,243. Their study suffers from some numerical 
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comp lications in the vicinity of trailina edge. 

Hummel [25] has made important rontnhution in revealing 
the complex flow structure behind a delta wing. He conducted 
experiments and presented pressure distribution, vortex structure 
etc of the flow around a delta wing of unit aspect ratio with an 
angle of attack of 20°. 

Many researchers have observed longitudinal vortices in 
various complex flow configurations. Naturally, the interaction of 
such vortices with boundary layer and its effect on heat transfer 

I 

was a subject of interest to them. Some examples are 
Tay 1 or-Gort 1 er vortices in boundary layers on concave curved 
surfaces, horseshoe vortices formed by an obstruction protruding 
from a surface and wingtip vortices impinging on a downstream 
surface. The embedded vortex is capabl® of strongly perturbing the 
boundary layer thickness and influpncing the heat transfer 
characteristics. In addition, longitudinal vortices usually 
maintain their coherence over a long streamwise distance. As a 
consequence, the heat transfer effects beliind a vortex generator 
are very persistent. Eibeck and Eaton [6], , Westphal and Mehta 
[26] worked extensively in this field with more focus to the 
turbulent layer. 

Eibeck and Eaton [6] conducted experiment on 
longitudinal vortices embedded in a turbulent boundary layer and 



resultant heat transfer effects. Longitudinal vortices are found 
to influence significantly the heat transfer behaviour. Local 
Stanton number increased as great as 24 percent resulting in a net 
increase in spanwise average heat transfer coefficient. Despite 
the presence of turbulent diffusion, the influence of longitudinal 
vortices on momentum and energy transport could be traced at a 
location as far downstream as 60 wing chords behind the delta 
wing. 

In the last decade researchers attempted to explore the 
aforesaid problems using numerical schemps. Rrockmeier et al C27] 
have computed the flow in a channel with huilt-in wing type vortex 
generator where full Navi er-Stokes equations were solved. 

Thomas, et al E283 have computpd low- speed laminar .flow 
over a low aspect ratio delta wing upto 40** angle of attack using 
an upwind biased finite volume algorithm. The differencing scheme 
used are second^order accurate and a multigrid algorithm is 
employed to promote convergence to steady state. The predicted 
results and lift coefficient of the wing f>ave remarkable agreement 
with experiments due to Hummel [2?] and others. At 40* angle of 
attack, a bubble-type vortex breakdown i«! evident in the 
computations. 

Biswas et al C5] have attempted to compute laminar 
mixed convection in a rectangular duct with built-in delta wing 



I using a modified version of MAC (Marker and Cell) algorithm. A 
number of computations have been performed with and without free 
convection for different values of Re. Hr ( Grashof number ) and 
angle of attack. Enhancement of heat transfer between the gas and 
the channel walls is evidenced. The augmentation is mainly due to 

4 

longitudinal vortices generated by the wing. Buoyancy induced 
secondary flow increases the vortex strength and improves the heat 
transfer still further. 

Westphal and Mehta C26] have studied the effect of 
oscillating vortex on a turbulent boundary layer. The meander is i 
simulated by forcing a periodic lateral translation of a 

half-delta wing vortex generator at a low frequency of 1 Hz and 
one half amplitude of 0.5 cm. The effect of the meander is found 

to flatten the vorticity contours at the stations where they are I 

I 

originally round. The Reynolds stresses are also found to be ' 
significantly affected. i 

The disturbance to the growth of thermal boundary layer 

i 

can also be accomplished by using a mul t i - louvered plate-fin ; 

surface. Investigations by Achaichia a Cowell [28] provides a 

detailed performance data for louvered fin surfaces. However, by i 

I 

using louvered fins, enhancement is obtained at the price of high i 
pressure drop. j 

The study of irreversibility in convective heat transfer j 


I 



13 


has become an important topic from thp standpoint of effective 
utilization of energy. It has been shnwn by Be)an ( C293, C503, 
[51] > that the irreversibility brought about by any heat transfer 
process results in penalty to useful power. In such analyses, the 
irreversibilities due to external interaction of energy and 
1 internal dissipative effects are s imu 1 taner^us ly taken care of. 
Estimation of irreversibilities of a process from the entropy 
production rate and the minimization of the same by a suitable 
adjustment of the governing input parameters is usually referred 
to as second-law or thermodynamic optimization of the process. 
Among the related works involving thermndynami c optimization, the 
analyses of heat transfer in cryogenic apparatus < Bejan and 
Smith, [52] >, study of counter-flow heat exchangers < Sarangi and 
, Chowdhury, [55] ) and swirling flow through a cylindrical duct < 
Mukherjee, Biswas and Nag, [54] ) are a few. The objective of the 
present article is to evaluate numerically the entropy generation 
of heat transfer processes in a rectangular duct involving 
augmentation by generating streamwise longitudinal vortices. Based 

i 

on an evaluation of irreversibility, this study draws some 
conclusions about the influence of type of vortex generators ( 
wings / winglets > on efficient use of energy. 



CHAPTER III 


STATEMENT OF THE PRf)PL FM 


Computation is performed in a chnnr^el which is formed by 
two neighboring fins < Fig,? ). An obstacle in the form of a 
delta-wing or wingiet pair < of zero thickness > is placed inside 
it. The base of the wing is fixed on the bof'tocn wall and the apex 
faces the incoming flow stream with an angle o+ attack. In the 
case of wingiet pair, one side is fixed nn the bottom wall and the 
trailing edge is free for each wingiet. Since symmetry prevails in 
the vertical central plane of the channel, the flow field in only 
half of the channel has to be computed. The dimensionless 
equations for continuity, momentum and energy may be expressed in 
the following conservative forms as 
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X c 
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<5> 

the 

above equations, velocities 

have 


nondimens i onal i zed with the average incnminq velocity U at the 

ctv 

channel inlet, all lengths with the channel height H. the pressure 


# 


/ 




Type of obstacle ( protrusion) 



Fig.5 Flow model for computation 
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3 

with pU and the nondimens ional temperature is defined as 0 = ( T 

av 

" T ) / ( T - T >. Boundary conditinns of interest in this 

00 W 00 

investigation are 

Top and bottom surfaces ; 

u = v = w = 0; T = T 

w 

Side wall ( z = B / 2 > and mid plane < z = 0 ) : 

w = ( dv/dz ) = ( du/&z ) = < 3T/9z > = 0 

At the channel inlet ; 

u = u<y) ,v=w=^0, T = T 

00 

At the exit, the boundary conditions for the dependent variables 
are obtained by setting the second derivatives of them in the flow 
direction equal to zero in order to ensure a smooth transition 
through the outflow boundary ; 

9 u _ 9 V _ 9 w _ 9 T _ Q 

* i” ~ a 2 ” ~ " ' 

dx 9x 9x ax 

No-slip boundary conditions fnr the velocities on the 

obstacles are used. The temperature of the obstacle is constant 

and equal to T . 

w 

While computing the case with stamping under the vortex 
generators, we have used periodic boundary conditions to take care 
of the holes due to punching on the flat plates. The period length 1 

i 

for these periodic boundary conditions is equal to the channel i 

t I 

height H C Fig.4<e> 1. However, the boundary conditions for the j 

fictitious cells on the bottom are : I 

I 

1 

I I 
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u. = u 

V = 0 . 5 ( V. , + V. , > 

v.l.k wJTRE.k t,l,k 

W. , = W. 

P = P 

v,l,k i,JRE,k 


and on the top : 


u. , = u. . 

i.ji>7XM,k v,2,k 

V. . = V. . 

, = 0 . 5 < V. , + V. , ) 

V,JRE,k \,,t,k v,JRE,k 

W. = W. . 

v,JXM,k t,2,k 

P = P 

i,JIM,k i,2,k 

The above two sets of boundary conditions are applied to 
all < i, k > locations on the hole. Such boundary conditions have 
also been discussed by Hirt, Nichols and Romero C??] 


Nusselt Number 

In order to have a quantitative distinction of the heat 
transfer performance the combined spanwise average Nusselt number 

H 




B [ q,4- q^ ] 


•a T^^<x,z)-T^(x)] dz + T^2 (x,z>-T^^<x)] dz K 

o'' o'' 


< 6 > 


'has been calculated at each longitudinal loca,tion to depict heat 
transfer at any axial position in the channel. The solution of 
energy equation provides with temperature value at each cell in 
the channel. The velocity quantities are known beforehand from'the 



solution of Nav 1 er-Stokes equations. These values are plugged in 
equation (6) to evaluate the combined 'soanwise average Nusselt 
number at each longitudinal section in th^ channel. 

The average Nusselt number for the entire channel is 
calculated as 


Nu 


e 



<7> 
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Entropy Generation 

The rate of entropy generation per unit volume is given 
by Kirkwood and Crawford C363 and Bird et al C37] as : ' 


S a L ( q . 7 T > (at 7 V > 

" T* T 


< 8 > 


The first term on the right-hand side of <8> may be 
written as ; ' 


. _i_ < , . 7 t > = r f?jf ] 

T® T* L Wx/ \dy/ \dzl J 


(9) 


The second term on the right-hand side of equation <8> 


is expanded as 
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After invoking incompressibility condition C equation (1>], 
equation (10) can be written as : 
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Use of equations <9> and (11) in equation (8) and 
nondimensional ization yields : 
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As expected, the irreversibility indicator Ns contains 

3 

two additive parts, one due to conduction in the presence of 
non-zero temperature gradient, and the other accounting for 


I 



viscous dissipation of mechanical pow^r throughout the flow. 

However, if the velocity and temperature fields are completely 

known, we can find out the values of nond imens i ona 1 entropy 

generation per unit volume < Ns > at each cell of the flow 

a 

I 

domain. After evaluating this quantity, we integrate Ns^ over the 
entire volume to get nondimens i ona 1 volumetric entropy generation 
in the channel as : 

I 

Ns NSg dX dY dZ <1?> 


Merit Function ' 

From this analysis, it is possible to evaluate the rate 
of energy transferred usefully as well as destruction of exergy 
due irreversibilities. 

If Q is the total rate of heat transfer, then 


a = h fT -f, ).2BL 

C > V o ■' 

or.Q*4NukLfT-f. } 
o V- V b 

In equation (14), a channel of aspect ratio 
considered, i.e., B / H = 2. 


(14) 
2 has been 


The equation (14) can be written as ; 



The rate of exergy transfer accompanying energy transfer 
at a rate of Q is given by Moran C'R] a- 


Q. = Q [ 1 - ] 

= a [ 1 - ^] 


(16) 


where T , the exergy reference environment temperature has been 

Ct 

considered as the ambient temperature T . The wall temperature T 

00 V 

has been considered as a suitable temperature at the surface where 

« 

heat transfer takes place. If S is the total rate of entropy 

■generation, the destruction of exergy is 
♦ 

I = T . S 

a 


or , I = 


^ T S 
T V 

V 


j] following the argument 3 


or , I 


. T S 

V 


(17) 


As such, the total rate of entropy generation, S. may be 


written as : 


/// S dx dy dz 

9 


or , S 


NSg d(HX) d(HY> d(HZ) 


(18) 


Invoking equation (1?) in equation (18), we obtain 


= H k Ns 


(19) 


Finally, equation (17) and equation (19) will yield 


I 


I 
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I = T H k Ns <20> 

r V 

A merit function is defined as the ratio of exergy 
transferred to the sum of exergy transferred and exergy destroyed. 


\ 


Q 


or , M 


a 

Q TT 


( 21 > 


a 

Substituting for Q from equation <15> and I from 

a 

equation <21), we obtain, 

4 < r-1 >® Nu < e - > 

M = S (22) 

4 ( r-1 )®'Nu <0 - 0, ) + Ns 

e V b L 

This merit function is now evaluated for various flow 
parameters in a channel using delta wing or winglet pairs as 

I 

vortex generators. The merit of the merit function lies in its 
simultaneous accountability of exergy transferred and its 
destruction caused by irreversibilities associated with exergy 
transport and momentum transport. Irreversibilities due to 
external interaction and internal dissipative effects are together 
taken care of by this parameter. However , it might be observed 
that equation (2i2> is another form of second-law efficiency < 


Moran , C?8] > . 



CHAPTER IV 


METHOD OF SOLUTION 

The nondimens i onal continuity, momentum and energy 
equations in conservative forms, written in terms of primitive 
variables C equations (1) - (?) ], are solved by using a modified 
version of Marker and Cell < MAC ) algorithm. The original version 
of MAC due to Harlow and Welch C?93 was modified by Hirt and Cook 
[40]. In the original MAC method, the pressure field was obtained 
by directly solving the Poisson's equation for pressure, whereas, i 
in modified MAC version, pressure values are calculated implicitly 

I 

from continuity equation by a pressure-velocity iteration process. ' 
A related technique developed by Chorin [41] involvec^ a 

simultaneous iteration on pressure and velocity components. 
Viecelli [42] has shown that the two methods as applied to the 

i 

MAC algorithm are equivalent. 

1 

The computational domain is divided into a number of 
rectangular cells of edge length 6Y and 6Z along the X, Y and 
Z directions respectively. Cells are denoted by an index < i,j,k > j 

implying the cell number as counted from the origin in X, Y and Z j 

j 

directions respectively. Staggered grid arrangement is used in 
which the velocity components are defined at the centre of the 
cell faces to which they are normal t Fig. 4(a> ]. The pressure 

and temperature are defined at the centre of the cell. 




Fig. 4(a) Three dirnentional staggered grid, showing 
the. locations of the discretized variables. 
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4 . 1 Computational Scheme 

MAC is a semi - impl i c i t scheme of solving complete 
unsteady Nav i er-Stokes equations where the advancement of velocity 
components with respect to time are obtained explicitly by 
calculating accelerations due to convection, diffusion and 
pressure gradient. After obtaining a provisional velocity field 
for a particular time step, the continuity equation is solved 
implicitly. Hence MAC is an implicit-explicit scheme. 

’ The complete Nav i er-Stokes equatibns are elliptic in 

space and parabolic in time. Hence their solution has to be time- 
marching. Again, since the equations are elliptic in space, we 

I 

need boundary conditions on all confining boundaries — even at 
the outlet. To start the computation, a guess field of velocity, 
pressure and temperature are assumed. To compute compressible 
flow, a certain .density field has also to be assumed ; and in. that 
case, we have to solve another equation, namely, the equation of 
* state. However, in our case, the flow-field is incompressible and 
we use constant density assumption. 

From the guessed velocity and pressure fields, the 

I 

corrected velocity and pressure fields are obtained by pressure- 
velocity iteration through continuity equation. Convergence of 
this iteration process ensures a divergence-free velocity field ! 
for the initial time-step. Now, this corrected pressure-velocity | 



field is used to calculate the velocii"y f i p ! d for the next time 

I I 

step by making use of Nav i er-Stokes en'j=itions. TPie advancement of 
the velocities for one time step <5r i <=; calculated by evaluating 
the accelerations caused by advectinn, diffusion and pressure 
gradient. The choice of the time increment value St is governed by 
stability criterion as discussed in spction 4,4. However, this 
explicit time advancement may not lead to a velocity field with 
zero mass divergence in each cell. In the subsequent second stage, 
adjustments are done to ensure mass mnservation in each cell. 
This is performed by manipulating the pressures as well as 
velocities in each cell through an iterative process, which, as 
mentioned earlier, is equivalent to solution of Poisson's equation 
for pressure. Details of pressure-velocity iteration procedure is 
discussed in section 4.5. The corrected pressure and velocity are 
now used to evaluate velocities for the next time step through the 
momentum equations. Thus the solution advances in time direction ‘ 
until the values of velocity components in between two consecutive 
time steps become identical implying steady state converged ' 
solution. ! 

The solution of the energy equation is, however, ; 
relatively simple. After the steady converged velocity field is , 
obtained, the velocity components are invoked in the energy j 
equation making it linear in nature. The solution of the i 

I 

i 

I 

i 


1 
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temperature field is then obtained by a Successive Over Relaxation 
< SOR ) technique. 


Discretization and Numerical Bou ndary Conditions 
The continuity equation is discretized by a simple first 
order differencing in such a way that the mass conservation is 
maintained. For a particular cell, mass divergence ( for 
incompressible flow > is given by : 


A. . = f U. - U 


) / ax 


f V. . - V. . , 1 / <5Y + r W - W ^ / 6Z 


The convective terms of Nav i er^-Stokes and energy 
equations are discretized by a hybrid scheme which combines upwind 
and central differencing to achieve the stability of upwind and 
formal accuracy of the central differencing < Spalding C45] )'. The 
discretization of one of the convective terms is shown below ; 


a(UV> 

dY 


fCV. ..+V. ,■) 

4 oY L v+l,j,k J k v,j,k \,j'±,k J 

+ af' V. ., + V. 1 f U. 

k v,j,k i.+l,j,k J k v,j,k '..j<i,k ■' 

- f V. . , + V. . 1 f LI. . , + u . , I 

k v,j-l,k t,+ l,j-l,k k v,j-l , k v,j,k J 

“ «f V. . ^ - V. . , ^ 1 f U. . , - U 1 ] 

k x,,j-l,k i-*-l,j-<-»k k v,j-j , k i.j.k ■' J 

Here a is a weighing factor. When « is zero, the above 



discretization will become a simple cet\tr?)l differencing and when 
a = 1, this is exactly the second upwind differencing. The choice 
of this factor ot is also dictated by stability criterion as 
discussed in 4.4. 

Discretization of temporal terms in equations <2) to <5> 

i 

is done by simple forward differencing. All the second order terms 
are discretized by central difference. For example, 


dU 


au _ /• , , , , n ^ 


ax* 


— - — f U. . - 2U + U 1 
<ax)* j 


An important aspect of solving differential equations is 
the application of proper boundary conditions. Correct 
implementation of boundary conditions in the solution code 
requires special treatment. For prescribing velocities on the 
boundairies of the computational domain, four kinds of velocity 
boundary conditions are needed, namely, rigid-wall or no-slip 
boundary condition, symmetric or free-slip boundary condition, 
inflow boundary condition and outflow boundary condition. At any 
boundary, one velocity component will be normal and the other two 
will be tangential to the boundary. The normal velocity is easy to 
prescribe when the boundary coincides with a cell edge. 


I 



At the inlet, the axial velocity component is given the 
desired input value or a well defined profile. The vertical and 
spanwise velocity components at the inlet are assigned zero value. 
For the treatment of the no- slip and froej-siip boundaries, we 
need to consider a layer of fictitious cells just outside the 
computational domain. On an no- slip wall, normal velocity 
components are directly set to zero, whereas, the tangential 
components on the wall are made zero by assigning the tangential 
velocities of the fictitious cells the same magnitude as those of 
the adjacent real cells but opposite sign. On a free-slip surface, 
the normal velocity is again zero, but here the tangential 
components in the fictitious cells have same magnitude and sign as 

I 

those for the cells just inside the computational domain. For 
outflow boundary there is no unique prescription, but the general 
idea is to choose the boundary condition m such a way that it 
has minimum upstream influence. For this purpose, a useful 
approach is to set the second derivatives of the velocities equal 
to zero. This seems to keep the flow going smoothly out of the 
boundary . 

For temperature boundary condition also, we have to make 
use of the imaginary cells. The temperature in the imaginary cells 
are set equal to the temperature values in t.he corresponding cells 
just inside the fluid region when a free-slip wall is concerned. 
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They are equal to ( 2T^ - T^^j^ > in car-e of a conducting wall 
and for the assumption of a linear variafion of temperature. 
However, a nonlinear temperature varir»tinn on the solid wall has 
been used in the present code, given by 

T1 = TW5 - F112 X Tl< l.J,K > ^ F?? x TK I,J,K > 
where, TW5 = 5 T 

V 

F112 =11/2 

F52 = ? / 2 

At the inlet, uniform tempeintuip profile has been 
employed. At the exit, the second derivative of the temperature in 
the axial direction is made equal to reto in order to ensure a 
smooth transition. 

There are special modules to take care of the boundary 
conditions on the wing generator. Subroutines BCD. BCOV and SCOT 
implement the velocity and temperature boundary conditions on the 
vortex generator. BCO is the subroutine where the geometrical 
parameters of the vortex generator are taken care of. From 
Fig.4<b) and 4<c>, the angle of attack of the vortex generator can 
be expressed as /3 = tan *< 6Y/6X >, aspect ratio A = 2b / t where 
b / 2 = < KC - KB >.i5Z and -C = < IB - lA >»6X / cos^. Fig.4<b> 
further reveals that the U and V components of velocity fall 
directly on the surface of the obstacle and, therefore, are set 
equal to zero. Implementation of no-slip condition for U component 


I 
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“ =■■ ' \]M 

Fig- 4(b) Relative location of velocity components 
■ and wing on x-y plane at Z=0 


I 



^c: 


velocity at the edge of the obstacle n»eds some manipulations 
which are incorporated in subroutine BCOV. From Fig.4fe>, it' is 


found that U 


U 


i*j»k 


ifj.k+l/'S 

U . 


0, or, in other words. 


The W component of boundary also needs manipulation for 
its zero value on the surface of the obstacle. From Fig.4<e), 

W. .. = - W. . , 

t#J»K v,j-l,k 




v,j,k 


- W, 


v-i, j , k 


W. . , = - W. .. 

V/j—A p k v,j,k 

w. , . . « - w. 

V“i,j , k v,j,k 

The temperature boundary conditions on the obstacle are taken 
care of in the subroutine BCOT and are diagfamat i cal ly shown in 
Fig.4<d>. The nondimens i onal temperature 9 has been denoted in the 
code by T. The obstacle boundary conditions for temperature are t 


T. 


V, j, k 


i. j. k 


= 2T - T. 

V v-l. J. k 


= 2T 


T. = 2T - T. . . 

i-l. j+i. k _ V v-1. J. k 

T. = 2T - T. , . ^ . etc. 

v-i. j+ i. k w v-8. j*i. k 


4. ? Pressure-Velocity Iteration 

Direct solution of Poisson's equation for pressure 
suffers from the inherent problem of prescribing pressure boundary 
























Fig. 4f Periodic boundary conditions for pressure and 
velocities on the top and the bottom wall 




conditions. In this respect, the modified MAC algorithm due to 
Hirt and Cook [403 has a considerable edge nver the original MAC 
method. The mathematically methodology nf the iteration process is 
described herein : 

If the velocity and pressure 'field are known at any time 
step, we can evaluate the velocity field for the next time step by 
explicitly calculating the acre 1 era 1 1 ons due to advection, 
diffusion and pressure gradients through Na v i er-Stokes eguations. 
As has been mentioned before, the explicitly advanced provisional 
velocities may not necessarily lead to a flow field with zero mass 
divergence in each cell, which means, at this stage the pressure 
distribution is not correct. Pressur»=> in each cell will be 
corrected in such a way that ther^ is no accumulation or 
annihilation of mass in the cell. The relationship between the 
explicitly advanced velocity component [ at fn+l>-th time step 3 
and the velocity at the previous time step [ i.e. n-th time step 3 
may be written as 


A# 

u 

v,j,k 


6tC P.^.. — ) 

_ y n ^ ^ J 

<5X 

+ 6t [ RFSIDU ] 


where f RESIDLI is the value of I - - 
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ax ■ 

- wr 


* Re 


at the ( i , j , k > cell evaluated with thp velocity values 


at the 
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n-th time step. ' 

On the other hand, the corrected velocity component ( 
unknown > will be related to the correct pressure ( also unknown ) 
in the following way : 

n.-! <5Tr ^ ) 

v.j.k i,j,k 

4. <5t r RFSIDLI Y" (Z4> 

From equations <23> and (24), one can write 


U 


n+l 
^ , j , k 


V, 4,k 






P. 


v + »,),k 


Si 6P' 


6X 


ax 


I i-k _ n 

where P. . , = P. . , - P. . , 

v,j,k t,j,k v,j,k 

, .. ar ap.'".. 

So. u.”'!’* = u.”’!'* + 

'.J.k v,j,k ^ 

Neither ap nor U are known explicitly at this 
i,j,k i,j,k 

stage so that one can be calculated with the help of the other. 
Calculations are done in an iterative cycle and we write / 


U. 


r>4-i 

^ f j ♦ ^ 


ar ap''.^ 

u, . , 

V ji J / K ^ 


In the similar way we can complete the following array 
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<5 t fiP 
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67 J 

The correction is done through 'continuity equation. 

Plugging in the relationship <25) into continuity equation, 
n-fi , , n+l . . n+l . . n+l . , ^ 1 

+ 


V — V W' — w 

i.,j,k i-l. j,k ^ v ,j,k i.,j-l.k ^ l,j,k i>j,k-l 






n+l 
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n+j 
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or. ap r - o> ( 
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(26) 
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1 


where is an over re 1 o -at i on factor which is 

introduced to accelerate the pressure correction process. Usually 


a value of = 1.7 is used. After calculating 6P'' , velocities i 


1 n 


each cell are corrected according to the equation set <25) and 
pressure in each cell is adjusted as 


, n+t 
wink 


4 6P^ 


(27) 


This process is continued until the velocity divergence 
in each cell vanishes. If the velocity boundary conditions are 
correct and a divergence-free converged velocity field is 
obtained, eventually correct pressure will be evolved in all the 
cells including the cells on the boundary. This feature of 
Modified MAC method has been discussed in more details by Peyret 
and Taylor C44]. However, it was also shown by Brandt, Dendy and 
Ruppel [453 that the aforesaid pressure-velocity iteration 
procedure is equivalent to solution of Poisson's equation for 
pressure . 


4 . 4 Numerical Stability 

The rectangular cell sides <5y , <5Y and SZ are related to 
the input parameters like angle of attack, wing aspect ratio, 
channel aspect ratio etc. Given the geometrical parameters of the 
flow model, «5X <5Y and 61 can be obtained, rjnce the mesh has been 
chosen, the choice of time incremen4 6t is dictated by stability 
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criteria which are governed by two restrictions. The first 
restriction is known as Courant-Fr i edr i chs-Lewy <CFL) condition 
which says that material can not move more than one cell in one 
time step. The difference equations assume fluxes only between 
adjacent cells. Mathematically, therefore, we should have 
<5 t < min [ <6X / juj), (6Y / jVj>, (6Z / |Wj) ] (28) 

where the minimum is with respect to every cell in the 

flow-domain. Typically 6r is choeter't equal to one-third to 
two-third of the minimum cell transient time. 

The second restriction is governed by 

gr i d-Four i er-number which states tha^- the momentum must not 
diffuse through more than one cell in one time step. A linear 


stability analysis shows that this limitation implies 


u 6r < 


6y^ <5z^ 


2 < ax*+ ay*+ <52^ > 

In nondimens ional terms this gives , 


1 6Z^' 

Z < <5X*+ «5Y*+ <57.® > 


. Re 


Finally, <5t is chosen such that it satisfies the above two 


inequal i t i es . 


It has already been pointed out that the term a in the 

discretization of convective terms of eqi.iations (1) ■ is an | 

! 

weighing factor'which gives the desired amount of upstream < donor 
cell > differencing. 0 = 0 gives space rotifred differencing and o ! 


I 
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* 1 gives full upstream 

provided the fluid is not 

one time-step. In general 

the maximum value of 

V 6t 


U 6t 

~ST' 


or 




or 


In other words 


or donor cell form, which is stable 
permitted to cross more than one cell in 
, a should hp rhnspn slightly larger than 

occurring ’in the entire domain. 


W 6 t 
61 


1 < o 


MAX 


r 1 u 

V 6x 

W <St 1 1 

L 1 -sx ’ 

"^Y 

-zr \ J 


An ot approximately 1.2 to 1 . times larger than the 
right hand number of last inequality is e good choice. Too large 
an a introduces unnecessary amount of rvimprical smoothening. < 
d i f f us i on- 1 i ke truncation error >. 


^ ^ S oftware Implementation 

The computations have been performed on a HP 9000 ( 850 

I 

I 

series > computer. The computer code consists of several modules. I 
*DELTAM' is the main program which uses a number of subroutines j 

* I 

each of which has a set of specific tasks to carry out. 

The given flowchart shows the overall structure of the 
code with its major communication links between different 

! 

! 

subroutines and main module. A short description of different j 

I 

indices and the main functions of different subroutines are given | 
below : 

IPARA This is an index which indicates whether the calculation | 



i 9 parabo lie < 


1 > or elliptic < 


0 > . 




ITYPE 

IRtlST 

INIT 

START 

RESTAR 

ICORRO 

JSTAM 

BCT 

BCTS7 

8CNS 

BCNSST 


This index denotes whether tlif' nbctaele is wing < = 1 > or 
wi tigl et < s Z > . 

This index denotes whether a re'^-tart of computing from an 
existing field variables are desired < = 1 ) or the 

computation should start from the very begining < = 0 ). 
This is a subroutine which initializes the entire 

calculations for a set of input vai'^iables such as Reynolds 
number, Prandtl number, grid sires etc. 

This subroutine creates the guessed field for the 

I 

dependent variables and initiates rnrnputat i on , 

This is a subroutine where the field variables are 

equated to a data field which was existing beforehand. 

This ipdex denotes whether the comDutation is for constant 
properties < = 1 > or vniiahln properties < = ? >. 

This is a parameter whicli denotes whether there is 

stamping on the plates ( - 1 ) or not < = 0 > 

This is the subroutine for temperature boundary conditions 
This is the subroutine for temperature boundary conditions 
when there is stamping on the plates. 

This subroutine provides boundary conditions for Navier- 
Stokes equations. 

This subroutine provides boundary conditions for Navier- 



CONTI 


VELALT 

TICOF^F^ 

OUTPUT 

ENERCiY 

STATE 

NSEQVP 


fiC|ua,tions when there t ‘^^^ti'p'nq on th<e plates. 

Ttiis is a main module winch solves continuity equation 
i mp 1 i c i t 1 y . 

This IS a subroutihe which storo's tfie current field 

variables as the solution for time r ~ r in order to 

o 

start a new .iteration for time r = t . 

o ‘ 

IN this subroutine, the value of time increment St and 
donor cell coefficient a are calculated from stability 
conditions in each cycle. 

This is a subroutine which prints tfie result files when 
the preassigned conditions for printing output are met. 
This subroutine solves the energy equation by an SOR or 
equivalent technique. 

This is the subroutine where the properties are updated 
through equation of state for variable property 
solut i on.NSEQCP In this subroutine velocities for the next 
time-step are explicitly evaluated through Nav i er-Stokes 
equations for constant property calculation. 

In this subroutine velocities for the next time-step are 
explicitly evaluated through Navier-Stokes equations for 
variable property calculation. 


I 


f 



4 - 4 . 


T I nr^An 


1hi«i 'tubrautine calculates tt>e time-gradients of the 
velocity components. 

In the following page, flow chart of the code is given. 
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results and PISCUSj^inM 


A large number of computations have been performed with 
delta-wing or a winglet pair as vortex generator under different 
operating conditions. Subsequently, a comparative performance of 
wing and winglet has been made through a second-law analysis of 
the heat transfer process in the channel with the wing or the 
winglet as the obstacle. Results have been obtained for steady 
incompressible laminar flow. Air has been taken as the working 
fluid; hence the Prandtl number of this study is 0.7. 
Incompressibility assumption is perfectly valid for flow of air in 
heat exchangers. Laminar flow assumption is also not for 
computational simplification. Generally, the passage between 
adjacent plate fins of heat exchangers are so small and the flow 
velocity is such that the flow is often laminar. 

Fig. 5 shows the cross-stream velocity vectors at 
different axial locations in the channel with delta wing as the 
vortex generator. The position of the chosen cross-sections with 
respect to the wing has clearly been shown. For slender delta-wing 

I 

in an infinite medium, leading edge vortices have a vortical 
structure and they diverge slightly. Now, when we envision a wing 
moving in an infinite medium, the wake grows longer which^ is 
basically a swirling flow supported by trailing vortices. However 
the flow in a channel with an attached delta wing shows no 
trailing edge vortices. Fig.? shows the generation of vortices and 












their elliptic deformation as they move •■’long the channel in 
presence of a delta-wing. Deformation tal-'es place due to the 
reduction in strength of the vortices V'liijch is brought about by 
the viscous resistance encountered by the vortical motion in 
course of its travel. 

Fig. 6 illustrates the static pressure distribution on 
two different cross planes in the channel with built-in delta 
wing. Here we have computed simultaneously developing flow with 
Reynolds number 1815. The nondimens ional static pressure ( p/pU^ > 

av 

contours depict that the vortical motinn i of free-vortex in 
nature. Our results follow the same qualitai-ive trends as those 
observed by Hummel C 26 ] in experiments. 

Fig. 7 shows the longitudinal velocity vectors along a 
vertical plane of the channel for the geometrical and flow 
parameters as those of Fig. 6. For this case of simultaneously 
developing flow, we have taken uniform volncity profile at the 
Inlet of the channel. The spiralling structure of the flow field 
due to the presence of the obstacle is evident from this figure. 

Fig. 8 shows the isotherms at different cross-planes in a 

I 

channel in presence of a delta winqlet pair. The vortex generators 
influence the temperature field strongly. The spreading of the 
isotherms over a wider region of the cross-section implies ^that 
the bulk temperature is higher in the downstream from the trailing 
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Re = 1815, Pr = 0.7 
3 = 20°, A = 1.25 
Of = 3.50 



(b) 


Fig. 6 'Distribution of static pressure at tvi/o 
cross- planes : (a) X = 3.5, (b) X = 5.63 
from the inlet of the channel 






Re = 1815,Pr = 0.7, A=1.25 
«=3.50, /3 =20° 



Fig. 7 Velocity vectors on a longitudinal section (x-y plane) at a distance 
away from the midplane of fhe channel. 
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of thR winglet pair. The irnportanrp of the secondary flow 

atui thr higf^er heat transfer is evident from this picture. 

As has been stated earlier, '^hall also consider the 

influence of punched holes beneath the wina generators on the heat 

tx'ansfer performance. Fig, 9 shows the comparison of cross-stream 

velocity vectors at different axial loration^ from the inlet of 

the channel for the cases without and with stemping on the solid 

walls. Due to the stamping on the walls a velocity field normal to 

/ 

the vortical motion is induced in the dr^wnstream direction. This 
induced downward normal velocity field redtjce®? the strength of the 
vortical flow and, as compared with the r:e.qe where there is no 
stamping, a decayed circulatory flow pattern is observed at ^ the 
same axial locations. 

Fig, 30 shows the heat transfer performance in a channel 
for a simultaneously developing flow in presnnce of a delta-wing. 
The figure also compares the heat tran^^fer for the cases of 
without and with stamping on the cfiannel walls. In the region of 
the wing ( from X = 2.92 to X - ^. 7 ^ the combined spanwise 

i 

average Nusselt number rises to a high value up to a region behind 
the middle of the wing and then takes a plunge. A small dead water 
zone exists in the immediate neighborhood behind the wing-wall 
junction which causes poor heat transfer at that location. Here, 
we are referring to the case with built-in delta-wmg in the 
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absencR of any stamping. However, in the downstream of the wing, 
heat transfer is increased remarkably as compared with the plane 
channel flow. As such, even for a very long channel < X = 10. Z8 >, 
the enhancement of the heat transfer at the exit of the channel 
is more than 46. 4Z percent. Enhancement in heat transfer is not so 
much when the effect of punched hole beneath the wing is taken 
into account. Due to the downward normal stream at the 
cross- planes ( as shown in Fig, 9 the strength of longitudinal 
vortices is reduced to a great extent and a spiralling flow with a 
relatively less vortex strength exists. The improvement in the 
heat transfer coefficient is relatively lower than that of the 
case of without any punched hole. However, at the exit of the 
channel, the improvement of heat transfer as compared with the 
plane channel is about 29.29 percent. This improvement is indeed 

' I 

better in all other axial locations behind the wing. 

Fig. 11 shows the influence of angle of attack of the 
winglet pair on combined spanwise average Nusselt number in ,the 
channel. At a nondimens ional distance of 4 from the inlet, for an 
angle of attack of 20’, we observe an enhancement of 50 percent in 
the combined spanwise average Nusselt number over the case of a 
channel without any obstacle. Now, at the same location, for an 
angle of attack of 26*. an improvement of about 6 percent in Nu^^ 

'over the case of the 20" angle of attack is obtained. However, at 
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thfi locatinn, an improvsrnent of l? pprc^nt over the cass of 
20" rtngle of attack is discerned for fi • 

F 1 g . 1 ?, compares the relativ* p*’ t- for mane e of delta-wing 
and winglet pair on augmentation of hf».nt francifar'. Heat transfer 
enhancement in the channel behind the vortey generators is clearly 
evident in both the cases. However, delt?’ wing produces streamwise 
vortices with higher strength wtnrt> brings about better 
improvement in heat transfer as compared with the case of the 
winglet pair. At a nondimensional distanre of 4 from the inlet, 
the enhancement in combined spanwise averaae Nusselt number Nu 

" oa 

for a built-in delta wing is more than '1.5 percent than that of a 
plane ctiannel. The enhancement for the bviilt--in winglet pair at 
the same location is about 21 percent. It may be mentioned that in 
case of built-in delta wing, a small dead-water zone exists in the 
immediate neighborhood behind the wing-body junction which 
eventually causes poor heat transfer at that location ( X = 2.5 ). 
The major difference between the wing and the winglet in this 
application is that the winglet has a free trailing edge whereas, 
the trailing edge of wing remains attached to the bottom plate. 
However, in the case of built-in winglet pair, due to the 
formation of trailing edge vortices, separation bubbles are not 
formed near the trailing edge and unlike built-in delta wing the 
zone of poor heat transfer is avoided. 
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1 




Fig, 12 Effect of type of obstacle on the 
distribution of cornbinecf spanwise average 
Nusselt nurober in tho nh'sjnnel 
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In order to estimate the flov*' losses in the channel with 
built in vortfty generators, the combinnd spanwise average friction 
coefficient has been defined as 

1 „dz 

■' y=H 

< 29 ) 

^ U* ( 2 B > 

Z av 
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r ^ 

[v 




dz + Z f ( 
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Fig. 15 shows the effect of varying vortex generators < 
winglet pair > angle of attack on combined spanwise average 
friction coefficient keeping the size constant. Increasing angle 
of attack has the effect of increasing vortex strength which 
increases resistance and consequently, higher value of combined 
spanwise average friction coefficient is obtained. At the exit of 
the channel, < x Re > for /? - 52° is 7.49 percent more than 
that for ft = 20° and < x Re ) for ft - 26° is percent more 

than that for ft = 20*. 

Fig. 14 compares the distribution of combined spanwise 

average friction cosfficient in a channel for built-in delta wing 

and winglet pair. As mentioned earlier, delta wing generates 

streamwise vortices with higher strength and as a consequence, ( 

C X Re ) in a channel is higher for delta wing. At the exit of 
t 

the channel, < x Re ) for a delta wing is about 146.6 percent 
more than that of the channel flow without any obstacle. In the 
case of built-in winglet pair, the increase in < x Re ? value 
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= 1.0, Hr = 0.39 

Without 1 

stamping 

^Without 

obstacle 


Location of the froiling edge 



Cf X Re 



Fig. 14 Effect of type nf nb‘=;tacle on the 
distribution of combined spanwise average 
friction coefficient in the channel 
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at the same location over the ( x Rr ) nf 3 plane channel is 
' nearly loo percent. It may be mentioned that the Reynolds number 
for this computation is 1500. 

1 ^ 19.15 shows the variation of merit function, M with 
Reynolds number ( in the range of 500 - ^000 ) for the built-in 
delta winglet pair. This figure shows tint tfm delta winglet pair 

i 

displays better performance than delta wing with regard to merit 

function. With increasing Reynolds number, the delta wing ; 

generates vortices with higher strenath which culminates in 

i 

steeper velocity gradients compared to the winglet pair. As a ■ 
result, frictional losses become high and consequently, the total 

entropy generation becomes much more pronounced for built-in delta ! 

, 1 

wing. For the case of heat transfer in a channel with built-in 

winglet pair, the increase in irreversibility associated with the ; 
increase in exergy transfer occurs at a relatively lower rate as 

compared to the case with a built-in delta wing. From the view 

1 

point of second-law efficiency, the use of built-in winglet pair | 
has edge over the built-in delta wing. Although with regard to I 
augmentation of heat transfer in the channel, the built-in delta 

\ 

wing shows better performance, it can be certainly argued that the ' 
built-in winglet pair is more effective with respect to energy 
transferred usefully. 






Because of non~avai labi 1 i ty of s'- oer imental data-base, 

the rnode 1 valiciation was not accomoli^hed extensively. However, 

one more comparison with respect to mean-Colburn-factor 

enhancement ( Aj ) using delta wing for two different angles of 

attack was performed. The experimental results were obtained from 

Kallweit (1986). In the experiment, the reference plate-fin area 

/ 

started from the location of ’the bes<=> of the wing; its 
longitudinal extent was 7,5 wing lengths and its lateral extent 
was 4 times the wing span. Identifying this zone on the bottom 
plate as the reterence area, experiments were conducted in a plane 
channel and in a channel with built-in rlelt-i wing at different 
angles of attack. Experimentally obtained mean-Colburn-factor 
enhancements < Aj ) over the referenrp aren due to the delta wing 
with angles of attack 20* and 50* are 1.1 x 10 ^ and 1.425 x 
10 ® respectively. Reynolds number for this experiment is 1815 and 
air is used as the flowing fluid < Pr-n.7 .>. Computations were 
performed in similar geometrical configurations by making use of 
our model < and hybrid scheme of discretization) and the 
corresponding mean-Colburn-factor enhancement values < Aj ) were 
found to be 1.0178 x lO'® and 1.7 x 10 ^ for the angles of 
attack of 20* and 50* respectively. The numerical model and the 
experiment corroborate each other reasonably. 
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CHAPTER VI 


CONCLUDING REMARKS ' 

Irreversibility, quantified by the rate of entropy 
generation , plays a significant role in convective heat transfer 
processes. The second-law analysis seeks to evaluate the rate of 
entropy generation in various processes and thereby, select the 
process which has minimum entropy generation rate and hence is, 
thermodynamically, the most viable onp. It has already been 
observed that significant improvements may be brought about in 
heat exchanger design by making use of wing-type vortex generators 
on the plate-fins. However, based on a comparative study of 
second-law efficiency, use of winglets appears to 
effective augmentation technique. 


a 


more 
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